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Abstract 

We investigate in detail the production of solar flare gamma ray emission above lOMeV via the 
interaction of high energy protons with the ambient solar atmosphere. We restrict our considerations to 
the broad-band gamma ray spectrum resulting from the decay of neutral pions produced in p-H reactions. 
Thick-target calculations are performed to determine the photon fluences. However, proton transport is 
not considered. Inferences about the form of the proton spectrum at 10-lOOMeV have already been drawn 
from de-excitation gamma ray lines. Our aim is to constrain the proton spectrum at higher energies. Thus 
the injected proton spectrum is assumed to have the form of a Bessel Function, characteristic of stochastic 
acceleration models, up to kinetic energies of a few hundred MeV, coupled to a Power-law m proton kinetic 
energy at higher energies. The detailed shape of the gamma ray spectra around 100 MeV is found to have a 
strong dependence on the spectral index of the power-law and on the turnover energy (from Bessel function 
to power-law). As would be expected the harder the proton spectrum the wider the 100 MeV featuie The 
photon spectra are to be compared with observations and used to place limits upon the number of partic es 
accelerated and to constrain acceleration models. 


Gamma ray production 

We require the production spectrum of 7-ray photons from the decay of secondary tt° produced in p-H 
collisions in solar flares. This is given by 
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where the integral over TL expresses the total number of pions per unit energy emitted by an average proton 
in time t, i.e. it allows for the proton spectrum to be modified by thick-target, losses as the protons traverse 
the flaring loop; j p (T p ) [protons/cm 2 /s/GeV] is the incident proton spectrum; A = 0»,,/m e )A (A 

is the Coulomb logarithm relevant to protons); S is the flare area, typical y cm , '' ■ 

are the total pion energy/momentum; T p is the proton kinetic energy and T P , min is the minimum kinetic 
energy required by a proton in the lab. system to produce a pion of kinetic energy 7V. The differential 
cross-section for the production of a pion with energy T p is given by 


<»> 

where (Ccr{T v )) denotes the inclusive cross-section for the reaction p + p -*■ *° + anything and (c/A/c/2,) 
is the normalised production spectrum of secondary . We have ignored contributions from collisions of 
protons with alpha particles and heavier ions. 

Pion production 

The production mechanism for tt° from p-p collisions is dependent upon the kinetic energy of the 
incident protons. Comparison of the various models with accelerator data at various energies show that 
the production mechanism can be split into three distinct regions (Dernier 1986): 
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T p < 3GeV : Isobar model (Stecker 1970) 

- pion production mediated by excitation of the A(1232) isobar which decays into a proton and a neutral 
pion. The distribution of pions from this decay mechanism is effectively isotropic in the lab. system. 

dN r jl/2 ~ m p 
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where 

w r is a normalisation factor. 

B is the isobar mass distribution given by a Breit-Wigner distribution, 
f is the normalised energy spectrum of pions in the laboratory system, 
s is the total energy availbale in the centre-of-mass. 

T p > IGeV : Scaling model (Stephens and Badhwar 1981) 

- pion production obtained via scaling arguments which imply asymptotic forms for the differential 
cross-sections 
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where E*(d 3 a* /dp 3 *) is the invariant cross-section for neutral pion production in p-p collisions (the asterisk 
denotes centre-of-mass quantities) and cosO max may be obtained from kinematic limitations. 


3 GeV <T P < IGeV : 

Between these two limiting energies neither of the above models fit the laboratory data very well and so 
we follow Dermer (1986) and extrapolate the normalised production spectrum linearly between these two 
limits. 


Incident proton spectrum 

Our choice of proton spectrum incident upon the target region is such that we have a Bessel function 
spectrum, characteristic of stochastic acceleration processes, up to some kinetic energy T p0 which is typically 
of the order of lOOMeV since stochastic acceleration is only valid in the non-relativistic regime and a power- 
law spectrum at higher energies, viz., 

jp(Tp) = A w A' 2 [2(3p p /m p aT) 1/2 ] T p < T p0 

= A N K 2 mZppo/m pQ T) l / 2 )(T p /T p o)- 6 T v > T p0 

where 

v4jv is a normalisation factor to be determined from comparison with observations . 

I< 2 is the modified Bessel function (cf. Ramaty and Murphy 1987). 
a is the acceleration efficiency. 

T is the particle escape time from the acceleration region . 

[typical value from fits to observations is aT = 0.02 - 0.04] 

T p o is the turnover energy from Bessel function to power-law . 
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Results 


In this section we display the form of our spectra for a variation of the spectral parameters, 6 (the 
power-law spectral index) and aT (the hardness of the Bessel function spectrum), for a turnover energy 
equal to the threshold energy for pion production, T p0 = 280A/eV, and the resulting photon spectra. We 
have also included the photon spectra obtained from a variation of T p o with aT - 0.04 and b - 4. 




Figure 1: Injected proton spectra, J p (protons/GeV/cm 2 /s) a) Bessel function with aT = 0.02 for 
T p < 280 MeV (pion production threshold) and power-law of spectral index 6 at higher energies, b) Bessel 
function spectra with range of spectral hardness aT (as a with turnover energy = oo). 
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Photon energy (GeV) 



in Figu^ri 2 a f BTsTeJ fun H PeC / fa ’ FP f ( photons / GeV / cm V«) resulting from the proton 
igures i . a) Bessel function/power-law spectra, b) Bessel function spectra. 


spectra show 






Figure 2c: Gamma-ray spectra 
those of Figure la with a fixed 6 = 


FpH (photons/GeV/cm 2 /s) resulting from proton spectra similar to 
-4 and aT = 0.02 and the turnover energy T p0 allowed to vary. 


Discussion and Conclusions 

Our theoretical calculations show that the detailed structure of the 70MeV feature depends strongly 
upon the form of the accelerated proton spectrum at energies from lOOMeV to a few GeV. The harder the 
proton spectra the wider the pion decay feature as one would expect. However, until a detailed comparison 
of our results are made with observation our conclusions are limited. There are a few observations which 
demonstrate a strong (possibly dominant) pion emission at energies of around lOOMeV. Events such as the 
SMM event of June 3 1982 (studied in detail by Murphy et al. 1987) and possibly a couple of events from 
1989 (Riecer - private communication) allow us some evidence for a neutral pion decay feature flattening 
the spectrum at ~70MeV. We intend to investigate these observations in more detail. Solar observations 
by the TASC spectrometer on The Compton Observatory should hopefully provide many more such data. 
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